Ground State Energy for Fermions in a ID Harmonic Trap 
with Delta Function Interaction 
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Conjectures are made for the ground state energy of a large spin 1/2 Fermion system trapped 
in a ID harmonic trap with delta function interaction. States with different spin J are separately 
studied. The Thomas-Fermi method is used as an effective test for the conjecture. 
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There is recent experimental and theoretical interest 
P, 0] in a model Hamiltonian for N one dimensional spin 
1/2 Fermions in a harmonic trap: 
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For N — 2, the eigenvalue problem of this Hamiltonian 
had been analytically solved [3|. We concentrate in the 
present paper on the ground state energy Ej for given 
total spin J as N — > oo. It is well know from group 
theory that for spin J the spin wave function is described 
by the two-row Young pattern [N — M, M] where M 
: Y — 2 J) 1/2 and the space wave function by its associated 
pattern \m. According to a theorem due to Lieb and 
Mattis [5j, when g is finite, 



Ej < Ei 



if J < J'. 



We notice the following facts: 
(A) As g — > +oo, the ground state energy Ej for any 
total spin J approaches [2| a limit 
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(B) For g = 0, the ground state wave function for spin 
J = N/2 — M can be taken to be a product of two de- 
terminants 0, [B| : 

* = det [uo(xi)ui(x 2 ) ■ ■ -un-m-\{xn-m)] 

x det [u (x N - M +i) ■ ■ .u M -i(x N )} ■ (3) 

where Ui(x) is the normalized eigenfunction of the har- 
monic oscillator. Its energy is 
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When J = Q, N = 2M and E = N 2 /A. Equations © 
and lHJ indicate that 



1/4 < Eq/N 2 < 1/2 for g > 0. 



(5) 



(C) For g -> -oo, when the total spin J = 0, M = N/2 
pairs are formed with spatial size of the order of —g~ x 
each. The internal energy of each pair is ~g 2 /A. Thus 
the M pairs contribute —Mg 2 /A to the total energy. Be- 
tween these pairs there are Fermionic repulsion as well as 
attractive delta function interaction. It is difficult to dis- 
entangle this complicated repulsion-attraction mix. But 
we observe that the Fermionic repulsion in absence of at- 
tractive g is given by J5J . So it is reasonable to surmise 
that as g — > — oo, the repulsion contributes X to Eq: 
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as g 



— oo, 



(6) 



is small corn- 
Equation ([6]) 



where N 2 /2 > X > 0. The X term in j6] 
pared with the other term as g — ► -co. 
holds for 1 < M < N/2. 

Equation ([5|) shows that for g > 0, En is of order N 2 
while ^ shows that Eq is of order M as g — > — oo. How 
can that be? This question leads to the following conjec- 
ture: 

Conjecture 1 : As N = 2M and N -> oo, the S /A 2 
g/\N curve approaches a limit: 



v ersus x 



15) and 



£o/A 2 - f (g/VN). 

show that 
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The schematic figure of conjecture is shown in Fig. 1 
and its generalization is as follows. 

Conjecture 2: For fixed J/N = b, the Ej/N 2 versus 



(11) 



x = g/vN curve approaches a limit 

Ej/N 2 -> f b (g/y/N), 

where 

/b(a;) — > 1/2 as x — > 
/ & = l/4 + 6 2 , 



+oo, 
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Simple Test of the Two Conjectures . It seems Table 1 The slope S[N - M, M] for different b. 

difficult to prove Conjectures 1 and 2. But there is a 
very simple but effective test: If Conjectures 1 and 2 are 
correct, then for large N and given J, the slope at g = 
is 

dEj_ = N3/2 d(E /N*) ^ Nm dm (15) 
dg d(g/VN) dt ' 



where t = g/y/N. Namely, the slope should be propor- 
tional to N 3 / 2 for large N. 

We calculate dEj /dg at g — by first order perturba- 
tion theory, 



dEj 
dg 



tpt(x)ij}](x)'tl)^(x)ip^(x)dx 



(16) 



where Vf an d ipi are fermionic annihilation operators 
and the bra &; ket dessignate the ground state at g = 0. 
According to ((3]) this ground state has N — M spin up 
particles and M spin down particles in single particle 
states uq, u±, .... Thus we have 



dEj 
dg 



N-M-l M-l 



X X] / Ui(x)u 2 j (x)dx. 



i=0 j=0 



Defining 



S[N ~ M, M] = N~ 3/2 



dEj 



dg 



(17) 



9=0 



we have 

S[N - M, M] = TV" 3 / 2 / p N - M (x)p M {x)dx, (18) 



where 



M-l 



Pm{x) = X u i( x )- 



(19) 
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The physical meaning of pm(x) is the density of the up- 
spin (or down-spin) particles at g = 0. 

The density Pm( x ) and the slope S[N — M, M] can be 
calculated in a straightforward way with a computer. We 
define a scaled density function: 



R M (y) = 



From ([19]) we have, as M approaches infinity, 

Rm{v) 
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> 1. 



(20) 



(21) 



Fig. 2 shows that Rm{v) tends to a limit as M — ► oo 
which wc shall calculate by the Thomas-Fermi method 
later in this paper. Here we list in Table 1 values of 
the slope S[N — M,M] for some given parameter b = 
J/N. One can see from the table that for each value of 
parameter b, S[N — M, M] rapidly converges to a limit 
as N increases, supporting both conjectures 1 and 2. 
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Thomas-Fermi Method For large values of M, we 
can evaluate pu by the Thomas-Fermi method (cf. Fig. 
3). When M is very large, the particles in dx have as 
their maximum kinetic energy (=maximum total energy 
minus the potential) ~ M — x 2 /2. Thus, the number of 
particles in dx is \/2M — x 2 dx/ir. I.e., 

as M — > oo. 



Pm{x) —> \JlM - x 2 /n, 
Thus from (f20"|) we have as M - 

Rm(v) -> V 1 - y 2 /v 



for 



\y 



< l. 



(22) 



(23) 



Considering (20-23), we obtain the slope S[N — m, M] 
as M -> oo from (JTSJ) 



lim S[iV-M,Ml 



1 2 2 



2VT + 26 
3^ 



(24) 



lim S[M,M] = ~ 0.135095. 



where 6 = J/JV, z = (1 - 26)/(l + 26), and F is the 
hypergeometric function. For J=0, reduces to, 



(25) 



The numerical values of (f24|) are listed in the last line for 
each value b in Table 1 for comparison. 
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FIG. 1: Ground state energy Ej vs. g curve for J = (Schematic). Equations (8) (9) and (10) refer respectively to regions A, 
B and C in this figure. 
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FIG. 2: Comparison of Rm(v) with its limit \/\ — y 2 /n. 
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FIG. 3: Thomas-Fermi Method. 



